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Abstract. It is shown that the one-dimensional or two-dimensional radially symmetric isothermal compressible 
Navier-Stokes system has no non-trivial global smooth solutions if the initial density is compactly supported. 
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This result is a generalization of Xin's work [6] to the isothermal case. 
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1 Introduction 



In this paper we consider the blow-up problem for the Navier-Stokes equations of compressible isothermal 
fluids 

f dtp + div(pu) = 0, 

I r)Jmi\ 4- HivYo?/ (9, iA — 11A11 — (\ 4- lAX7r\ivn 4- X7 P( n\ = 



d t (pu) + div(pu ® u) - p,Au - (A + ^i)Vdivu + VP(p) = 0, 
with initial data 

p(x,0) =po(x), u(x,0) =u (x), l<n<2. (1.2) 

This system describes the motion of viscid gas. u G K™, p, p and A denote the velocity, the density, the 



■ coemcient of viscosity and the second coefficient of viscosity, respectively. The pressure P{p) satisfies 

P(p) = ap (1.3) 

for a positive constant a. 

Mathematically, system is a quasilinear transport-parabolic system. There are many important 
;_h ' questions regarding this system. One of them is whether the solutions will develop singularities from 

smooth initial data. To our best knowledge, Xin [6] first proved that when the initial densities are com- 
pactly supported, any smooth solutions to the compressible Navier-Stokes equations for non-barotropic 
flows in the absence of heat conduction will blow up in finite time for any spatial dimension, and this 
feature also holds for the isentropic flows (i.e. P(p) = ap 1 with 7 > 1) in one dimensional case. Cho 
& Jin [T] extended Xin's work [B] to the case of fluids with positive heat conduction. Recently, under 
an additional assumption that one of the components of initial momentum is not zero, Tan & Wang 
[5] gave a much simpler proof of the result of Cho & Jin pQ. If the initial data do not have compact 
support but rapidly decrease, for n > 3 and 7 > Rozanova g] proved that any smooth solutions to 
the compressible Navier-Stokes equations for the non-barotropic flows with positive heat conduction still 
blow up in finite time. 

It is worth mentioning that very recently, Huang, Li, Luo & Xin [3] proved that any smooth radially 
symmetric solutions to the two-dimensional isentropic fluids with compactly supported initial densities 
blow up. And they also proved in j3 a an interesting result that, under some compatibility condition for the 
initial data, the small smooth solutions to the three-dimensional isentropic compressible Navier-Stokes 
equations with compactly supported initial densities are globally well-posed. 

Besides the non-barotropic ones and the isentropic ones, the isothermal flows are also of great physical 
significance. From the view point of mathematical structures, the isothermal case could be looked as the 
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endpoint of the isentropic case. Roughly speaking, we show that if the initial data are nontrivial and 
radially symmetric, and the initial densities are compactly supported, then the smooth solutions to 
will blow up in finite time. More precisely, we prove that: 

Theorem 1.1 Assume that 

2 

/i > 0. A+-/i>0. (1.4) 

n 

and that po{x) has compact support. Suppose that one of the following conditions holds true: 
(i) n = I: 

(ii) n = 2, and the initial data are spherically symmetric, i.e. p = p(\x\, t), u = u(\x\,t)x/\x\. 

Then the solution (p,u) € C 1 ([0, T], H m (W n )) (m > 2) to system with nontrivial initial density will 
blow up in finite time. 

Remark 1.1 One key assumption in Theorem \l.l\ is that the initial density is compactly supported. This 
assumption could be relaxed to fast decay at spatial infinity like the paper CT. But it seems pretty hard to 
get some kind of blow-up result for the initial density with positive lower bound. 

Remark 1.2 Our method also works for the isentropic case and non-barotropic case. The proofs essen- 
tially are the same. One reason is that all these systems have the similar energy conservation and the 
same momentum equation. 

To prove Theorem ll.il we use the so-called function method which is typical in the theory of blowup. 
The idea is to integrate the equation suitably and then try to deduce blowup. For the current theorem, 
the proof roughly goes as follows. First we show that the density keeps being compactly supported, 
which is similar to that of [6]. Then we integrate the momentum equation by the weight x, which gives 
an integral identity. It consists of three terms. One important observation is that, as time grows, two 
terms keep bounded, one term grows linearly. This fact implies no nontrivial global smooth solution 
exists. The details will be presented in Section 2. 



2 Proof of Theorem 11.11 

In this section we give the proof of Theorem ll.il which consists of four steps. 
Step 1. The density p is compactly supported all the time. 

Because the initial density po(x) has compact support, there exists a constant R > such that 

supppoO) Q B R . (2.1) 
Denote by x(t;x) the particle path starting from x, i.e. x(t;x) satisfies the following equation 

dx , 

— = u\x, t), , 
dt V ' h (2.2) 

x(0) = x, 

where u is the velocity. Denote by f2 t the closed region that is the image of Br under the flow map (|2.2[) , 

ilt := {(a;,i)|:r = x(t;x),x £ Br} . 

Note that on the particle path x = x(t\x), the density p satisfies a homogeneous ordinary differential 
equation. Thus, if p is zero at some initial position x, then on the particle path starting from x, p will 
be zero all the time. By (|2.ip . we see that 



p = o in ni 

Subsequently, using the momentum equation of (II. ip . we get 

p,Au + (A + ^)Vdivu = in f2jr. (2.3) 
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In the one-dimensional case, owing to (|1.4[) . we deduce from (|2.3p that 

U xx = in Q£ . 

This in combination with the condition u(-,t) £ H m (M}) implies that 

u(x, t) = in J7j . 

Therefore, in view of (|2.2I) . we derive 



= C Br. 

Lai case, because m(x, i 

have from (I2.3[) that 



x 

In the two-dimensional case, because u(x,t) = - — -u(\x\,t) for some radially symmetric function u, we 

\x\ 



a , , ( - J =0 in <»; . 



(7), 

Here r := x|. Using the condition u(-,t) £ H m (M. 2 ) again, it follows that 

u 

u r H — =0 in fl 1 : . 
r 

One can easily compute that the general solution to this equation is 

/ \ C(t) 
u(r,t) = 

r 

where C(t) is a constant that only depends on t. Because u £ C([0, T], i 2 (R 2 )), we get for a.e. t £ [0,T], 

+00 

\u{r,t)\ 2 rdrdS < +00, 



U12 Jo 
which gives 

u(r, t) = in fi£. (2.4) 

Therefore, fit = ^ 

S'iep 2. Averaging the equation by a weight. 

We multiply the momentum equation of ([l.ip by weight x and integrate the resulting equation on the 
whole space to derive 



pu ■ xdx 



PoUq ■ xdx = I I (p\u\ 2 + nap) dxdr, 
Jo Js," 

where we have used the isothermal condition (|1.3[) . In view of the fact that 

p(x,t)dx = / po(x)dx := itlq, (2-5) 



and that supp x p(x.t) C Br, we then get 



pu ■ xdx — j PqUq ■ xdx = / / p\u\ dxdr + namot. (2-6) 



Obviously, the right-hand side of (12.61) grows linearly. Next we show the left-hand side of (|2.6I) is bounded. 
This will give the contradiction with the induction hypothesis that the Navier-Stokes system has a global 
solution. To get the boundness of the left-hand side of (|2.6|) . we need to prove energy conservation. 
Step 3. Energy conservation. 

Multiplying the density equation of by ln(p + e), where e > is a small constant, we get 

Pt ln(p + e)dx = — div(pit) ln(p + e)dx. 



It is easy to see that 

p t ln(p + e) = [pln(p + e)] t - p[hi(p + c)]t = [pHp + e)]t - ^ = IpHp + c)]t - Pt + 4Hp + e)]t, 



3 



and that 

div(pu) ln(/7 + e)dx — 1 ' 



p + e 

j u ■ \7pdx + u ■ V[eln(/9 + e)]dx 
/ u ■ Vpdx — / divu ■ e ln(p + e)dx 

J Br J Br 



[pln(p + e) + eln(p + e)]dx = / u ■ \7pdx + / divu • e ln(p + e)dx. 



Owing to (|2.5p . we then have 
~dt 

Let e -> 0, note that p G C([0, T]; H m (R n )) for 1 < n < 2, the function p ->• plnp is continuous with 
respect to p G [0, +oo) and that |eln(p + e)| < C|elne| — > 0, it is easy to calculate that 

/ p\npdx= / / u-\7pdxdr+ / p^hipodx. (2-7) 

J Br Jo J Br J Br 

Subsequently, multiplying the momentum equation of (jl.lj) by u and using the equation of density, we 
obtain 

— / p\u\ 2 dx + p \Vu\ 2 dx + (A + p) I \divu\ 2 dx + a u-Vpdx = Q. 
2dt J R n J Rn J Rn J Rn 

Subtracting this identity from (|2.7p gives that 

— f p\u\ 2 dx + a j p\npdx + p f I \\7u\ 2 dxdr + (A + p) f j \divu\ 2 dxdr 

2 JBr J Br JO J Br JO J Br 

= o / Po\uo\ 2 dx + a / p Q \np a dx. 

1 JBr JBr 

Owing to the right-continuity of plnp at p — 0, we get 

i-T f 

\S/u\ 2 dxdr < C. (2.8) 



'0 JBr 

In the one-dimensional case, by Newton-Leibniz' formula, we get 

\\u-xf Loa <R 2 \\u\\ 2 L ^ <R 3 / \u x \ 2 dx. 

J-R 

X 

In the two-dimensional case, because u(x 1 t) — —u(r,t), it is easy to calculate that V^itl 2 — \u\ 2 /r 2 + 
\d r u\ 2 . Then using the Newton-Leibniz' formula and (12.41) . 



\\u{x,t)-x\\l T =\\ru{r,t)\\l ? <2R^J \u(r,t)\ 2 dr + j r 2 \d r u(r, t)\ 2 dr} 

'Xf R ^T^-rdr+ f R \d r u(r,t)\ 2 rdr) 



< 2R 2 < 



R 2 { 

< 7T" / \V x u{x,t)\ 2 dx. 
2lT Jb r 

Consequently, 

\pu-x\dxj < m 2 \\u ■ x\\ 2 L ^ < C \\7u(x,t)\ 2 dx. (2.9) 



IB R 

Step 4- Derivation of the blow-up 
By ([ZI]), $2JB and (pT8]l . we have 

r T , r , r J 

a 2 m 2 r,T 3 < 2 



^ J pu ■ xdx^j dt + 2 j ^ j PoUq ■ xdx^j dt 
IT. 

Therefore, it follows that the smooth solutions to (jl.lj) will blow up in finite time. □ 
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